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ABSTRACT 

o 

o _ 

I Galaxy groups likely to be virialized are identified within the CN0C2 intermedi- 

' ate redshift galaxy survey using an iterative method. The resulting groups have a 

OJj' 1 

^ ' median velocity dispersion of about 200 kms . The virial mass-to-light ratios, us- 

■ ing k-corrected and evolution-compensated luminosities, have medians in the range of 
150 — 25O/1M0/L0. The number- velocity dispersion relation is in agreement with the 
low-mass extrapolation of the cluster normalized Press-Schechter function. The two- 
point group-group correlation function has rg = 6.8 it 0.3/i~^ Mpc, which is larger than 
the correlations of individual galaxies at the level predicted from n-body calibrated halo 

(•~^ I clustering. We conclude that the global statistics of groups are in approximate accord 

' with dark matter halo predictions. The groups are stacked in velocity and position 

■ to create a sample large enough for measurement of a density and velocity dispersion 

■ profile. The resulting stacked group contains about 1000 members above a well de- 
fined background distribution. The stacked mean galaxy density profile falls nearly as 

Q^i a power law with r~^'^ and has no well-defined core. The projected velocity dispersion 

is examined for a variety of samples with different methods and found to be either flat 
^ I or slowly rising outwards. The combination of a steeper-than-isothermal density profile 

^ I and the outward rising velocity dispersion implies that the mass-to-light ratio of groups 

rises with radius. The M/L can be kept nearly constant if the galaxy orbits are nearly 
circular, although such strong tangential anisotropy is not supported by other evidence. 
The segregation of mass and light is not dependent on galaxy luminosity but is far more 
prominent in the red galaxies than the blue. The M/L gradient could arise from orbital 
"sloshing" of the galaxies in the group halos, dynamical friction acting on the galaxies 
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in a background of "classical" collisionless dark matter, or, more speculatively, the dark 
matter may have a true core. 

Subject headings: cosmology: large-scale structure, galaxies: evolution 

1. Introduction 

Small groups of galaxies are important cosmological indicators of the distribution and proper- 
ties of the dark matter in the universe. They occupy the mass and velocity range between individual 
galactic halos and the large halos of rich galaxy clusters (Abell 1958; Burbidge & Burbidge 1960; 
Gott & Turner 1976; Hickson 1982; Ramella Geller & Huchra 1990; Nolthenius & White 1987). 
The RMS velocity dispersion of groups is only somewhat larger than that of individual galax- 
ies, but groups have the advantage that visible galaxy tracers extend throughout the dark matter 
halo. At intermediate redshift, groups are suitable targets for X-ray observation and weak gravita- 
tional lensing studies which are complementary probes of their contents.. Consequently groups can 
be used to probe the properties of the dark matter on scales and at velocities much smaller than 
can be examined in galaxy clusters. 

The theory of structure growth in the universe is based on the paradigm that the dark matter 
consists of collisionless particles that only interact via the gravitational force. Cold Dark Matter 
is a specific form of this hypothesis that has been subjected to intensive theoretical study. A 
particular strength is that the properties of virialized halos can be predicted from a given density 
perturbation spectrum to full nonlinearity via simulations and various analytic approximations. 
These predictions have been tested with varying degrees of success against the dark matter halos of 
individual galaxies and rich clusters, but are less examined on intermediate scales. The intermediate 
scales are interesting because they are at much higher phase space densities than massive galaxy 
clusters, yet their central dark matter densities are not overwhelmed and altered by the baryonic 
matter, as is the case in for most normal galaxies. 

The CDM theory gives specific predictions of the statistical properties of the dark halo 
population and the mean internal properties of individual halos. The Press-Schechter theory 
(1974, hereafter PS) predicts the numbers of halos as a function of mass or velocity disper- 
sion. At low redshift, suitably selected groups have a population volume density in accord 
with the cluster- normalized PS prediction (Moore Prenk & White 1993; Girardi &; Giuricin 2000). 
A second global statistic is the clustering of dark matter halos which is predicted us- 
ing analytic approximations which have been compared to n-body results (Mo &i White 1996; 
Jing 1998). This biased "peaks" theory predicts a slow increase of clustering strength with 
halo mass. The internal density structure of the halos is found in simulations to have a 
power law cusp, r~^'^ to r~^'^ which asymptotically steepens to approximately beyond 
the virial radius (Dubinski &; Carlberg 1991; Navarro, Prenk & White 1996; Moore et al. 1999a; 
Avila- Reese, Firmani, Klypin Sz Kravtsov 1999). The goals of this paper are to find a collection of 
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virialized groups and then to compare the predictions of the global statistics and internal structure 
of intermediate mass dark matter halos to the observationally derived properties. 

The paper is organized as follows. The next section describes our approach to identifying groups 
in a redshift catalogue. Section 3 gives an overview of the Canadian Network for Observational 
Cosmology's field galaxy redshift survey (CN0C2) and the virialized groups that we find. In 
§4 the number-velocity dispersion relation and the two-point correlation function are computed 
and compared to dark matter halo predictions. In §5 we "stack" the groups on their centers 
and measure the mean projected density distribution and projected velocity distribution. In §6 
we model these as projections of simple 3D systems to derive the mean mass density profile and 
mass-to-light ratio as a function of radius which leads to the discovery that groups have a rising 
mass-to-light ratio. We examine the rising mass-to-light ratio for various galaxy sub-populations to 
search for systematic trends that might point to its physical origin. We conclude with a discussion 
of the possible implications of these results and a short set of empirical conclusions. We use 
Hq = lOO/ikms"^ Mpc~^ throughout this paper and adopt JIm = 0.2, = as our reference 
cosmological model. The distances and transverse lengths would be about 8% larger in an CIm = 
0.3, A = 0.7 cosmology at the median redshift. 



2. Finding Groups in Redshift Space 

A group is defined here as a collection of three or more galaxies, above some minimum lumi- 
nosity, that meets a set of positional requirements designed to minimize chance associations. One 
might seek, say, bound groups, such as the low overdensity Local Group, or, virialized groups, 
which are collapsed and hence quite dense. Virialized groups are a high density subset of bound 
groups and are the aggregates of interest in this paper. Unfortunately, in redshift space there is 
a fundamental degeneracy between position and line of sight velocity. Consequently, the precise 
galaxy membership of a group found in redshift space is always a statistical issue. 

Most group search methods are based on the friends-of-fr lends (f-o-f) algorithm used by Huchra 
& Geller (1982). This is an important method that gives unique groups, independent of starting 
galaxy. The f-o-f algorithm can be tuned to yield groups of varying overdensity. The algorithm 
starts with any galaxy as the beginning of a trial group. All galaxies closer than some maximum 
distance (discussed below) are added to the group. Then, each of the new group members is in turn 
used as a center to search for its neighbors to add to the growing group, continuing until no more 
new neighbors are found. Then one proceeds to a previously un-examined galaxy to try to start a 
new group. This process continues until all galaxies have been examined for neighbors. Groups of 
one and two are then deleted from the catalogue. 

The f-o-f algorithm has two parameters in redshift space: a maximum separation in projected 
radius, r™"^, and either a maximum separation velocity, Av™"^, or co-moving redshift space dis- 
tance difference, r™"^, required to join the group. The Au™"^ and r-^"^ parameters are related 
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through Av = H{z)rz/{1 + z). We use Av for kinematic measurements and for group finding. 

The r™"^ and r^°^ parameters need to be mutually adjusted to take into account the mean volume 
density of the survey, n{z), so as to produce an overdensity with respect to the field chosen on the 
basis of an experimental goal with an allowance for redshift space blurring. The resulting galaxy 
overdensity in the cylindrical redshift space search volume is, 



no 27r[r^"^]2 rf n(z) ■ ^ ^ 

The r^"^ parameter needs to be chosen in relation to r^"^ with some care. One approach is to 
tune it to n-body simulation results (Nolthenius & White 1987). However, an ever-present issue is 
that the resulting group sizes and velocity dispersions tend to correlate with the r^"^ and r^"^ 
parameters. 

Here we are interested exclusively in virialized groups, suggesting we devise a variant of the 
basic f-o-f algorithm which selects groups that appear to be sufficiently dense that they will quickly 
virialize. In configuration space virialization demands a mean interior density of approximately 
178f2°-^^/9c (Eke, Navarro &; Prenk 1998), equivalent to nearly 350po in a low density, flat cosmology. 
A conventional approximation is that such groups will be contained inside the radius r2oo which can 
be estimated from the virial theorem as r2oo = y/3o'i/[WH{z)], where ai is the line-of-sight velocity 
dispersion and H(z) is the Hubble constant at the redshift of interest (Carlberg et al. 1996). For a 
given o"! the virialized group members will be, on the average, contained with approximately 1.5r2oo 
if the mean density is falling like r~^. This relation between velocity and projected separation 
immediately suggests a natural range for the parameters. The field pairwise velocity dispersion is 
approximately 300 kms~^ (Davis & Peebles 1983), equivalent to a single galaxy random velocity of 
about 200 kms^^. This velocity dispersion is generated in groups which suggests that an average 
^200 will be about 0.3h~^ Mpc at z ~ 0.4. 

The f-o-f algorithm provides a set of trial groups whose properties are fixed by the input 
link distance parameters such that many of them may not be virialized. Virialized groups have a 
minimum overdensity of about 200pc- Hence, for each trial group we estimate a velocity dispersion 
which is then used to calculate r2oo- Galaxies beyond a distance to the group center of 1.5r2oo 
are discarded. The remaining galaxies are used to recalculate the velocity dispersion. This can be 
iterated until the group converges to a stable set. Some trial groups quickly drop to only one or 
two members and hence no longer qualify for group status. On the other hand, if we choose a very 
large starting value of r^"^ or Au"*"^ a few groups will percolate over very large structures. A 
minor complication is that we must identify a group center. 

The details of our group finding algorithm follow. (1) Pick a cosmology for the analysis 
{Hq = 100 kms^"*^ Mpc^^, = 0.2, IIa = 0). (2) Set the sample's redshift and absolute luminosity 
limits (k-corrected and evolution-compensated at a mean rate of one magnitude per unit redshift) 
of = —18.5 mag, no initial redshift limits) which defines a galaxy sample for all further 

operations. (3) Pick an r^^^ (our standard groups use 0.25^~^ Mpc) and r^"^ (5 h^^ Mpc for our 
standard groups). Center a cylinder of radius r^"'^ and forward and backward extent of r^°^ on 
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each sample galaxy and count the number of sample galaxies. To create a background estimate 
we randomly draw points from an n{z) fitted to the full sample and count the number within the 
sample cylinder. If the search radii initially give less than three neighbors, then multiply smoothing 
lengths by 1.5 and repeat for this sample galaxy. At this stage a new subsample is defined requiring 
that the local overdensity relative to the smooth sample have some specified minimum value. In 
this paper we only impose the requirement that the local overdensity be positive. (4) Select the 
highest density ungrouped galaxy and begin to find a new group. (5) Standing on each new group 
member in turn, add to the group any galaxy in the minimum overdensity subsample that is closer 
than r™"^ and rj"*^. Repeat this step until no new galaxies are added. (6) This f-o-f group defines 
the starting group for the virialized group iteration. (7) For the trial virialized group, determine the 
geometric selection function weighted mean x, y, z and ai. Trim galaxies beyond, or add galaxies 
from the f-o-f list within, Vp = 1.5r2oo and Av = Sai. Repeat step 7 four more times, requiring 
that the last two iterations have an identical result. (8) Drop single galaxies and pairs from the 
catalogue. 



3. Groups in the CNOC2 Survey 

The Canadian Network for Observational Cosmology Field Galaxy Redshift Survey (CN0C2) 
was undertaken primarily to study the dynamics of galaxy clustering at intermediate redshift. The 
survey methods and catalogues are fully described in Yee et al. (2000). The survey covers a total 
of about 1.5 square degrees in four patches spread around the sky for observing efficiency and to 
control cosmic variance. Galaxy redshifts are obtained over the redshift range to 0.7, with the 
unbiased spectroscopic sample extending between redshifts 0.1 and 0.55. The catalogues contain 
approximately 6000 galaxy redshifts with an accuracy between 70 and lOOkms^^, along with 
UBVRI photometry. The groups are constructed from these catalogued galaxies. The luminosity 
and clustering evolution of the sample as a whole has been previously discussed (Lin et al. 1999; 
Carlberg et al 2000b). 

The CN0C2 sample has mji < 21.5 mag and well defined spectroscopic completeness weights. 
The average redshift completeness is about 45%, with nearly 100% completeness 3 magnitudes 
above the limit and about 20% at the limit. Within our primary redshift range, 0.10 to 0.55, the 
redshift completeness to the fiux limit is higher than for the sample as a whole. From our luminosity 
functions (Lin et al. 1999) we calculate the probability that we will successfully obtain a redshift 
as about 75% with the other 25% being for galaxies at redshifts out of our primary redshift range 
(Yee et al. 2000). Together these imply that about 60% of the galaxies above the flux limit in the 
redshift range 0.1 to 0.55 have measured redshifts. If a group contains 3 galaxies within the survey 
limits, then the probability from the cumulative binomial distribution that we will obtain 3 redshifts 
is (0.6)^ or 0.216. If the group contains 4 eligible galaxies then the probability we will obtain 3 
or more redshifts rises to 0.47. The probabilities of obtaining 3 or more redshifts in groups of 5, 
6, 7 and 8 eligible members are 0.68, 0.81, 0.90 and 0.95, respectively, assuming no complications 
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from geometric selection (Yee et al. 2000). Our average group contains 3.8 galaxies with redshifts. 
We conclude that our roughly one-in-two sampling allows us to detect about half of the three or 
more members groups that are present. This level of completeness has no bearing on most of our 
analysis so we normally do not attempt to compensate for this effect. 

There is a small effect due to groups at the boundaries of the surveyed region. Examining a 
typical field we find that no more than 20% of the groups have any of their r2oo area beyond the 
boundary and no more than half of that. Consequently approximately 10% of the group members 
are "missing", which will have little effect on the velocity dispersions, and cause roughly a 10% 
diminution of the group masses, through the virial radius, and a similar effect on the total light. 
In the presence of much larger random fluctuations due to the small numbers in the groups this is 
not a major concern. 

The only two parameters that turn out to have much of an impact on the groups are r^°^ 
and r""*^. Our density requirement forces the "raw" groups to have a velocity dispersion strongly 
correlated with the search distance in the redshift direction, 7-™"^. The effect of the iteration on 
the group velocity dispersion is shown in Figure 1. The initial velocity dispersion is calculated for 
the f-o-f groups and is hence very strongly correlated with the chosen r^°^. The iteration to find 
virialized groups allows the ai values to relax to more appropriate values. However some less than 
ideal groups do persist which will lead us to continually consider alternate samples throughout this 
paper. Most of the groups converge to a stable membership in one or two iterations, with only 
1-2% discarded due to failure to converge in four iterations. 

Figure 2 shows the redshift distribution of the groups for four group catalogues with increasing 
j.max_ Beyond redshift 0.45 the sample becomes incomplete as the flux limit of niR = 21.5 mag 
causes galaxies to slip below the sample absolute magnitude limit = —18.5 mag. Both this 
figure and a detailed examination of the lists of groups shows that the sets of groups have a 
large overlap, independent of the search parameters. That is, this indirectly indicates that the set 
of group centers is relatively insensitive to the group finding procedure. In the next section we 
examine a variety of statistics to form a basis to select some group catalogues as best suited to 
various analyses. 



4. Global Properties of Groups and Halos 

We will use three global properties of the groups to assess the degree of correspondence of var- 
ious group catalogues to dark halos. These are the mass-to-light ratio distribution, the abundance 
as a function of velocity dispersion and the clustering properties. The virial mass to R band light 
ratio, Mvt/Lr, of groups is an indicator of the value of CIm- The number density of dark matter 
halos as a function of their one-dimensional RMS velocity dispersion, n{ai), is an important test of 
the CDM clustering spectrum but here is used as a guide to whether the number of high velocity 
dispersion groups is reasonable. The group-group auto-correlation as a function of mean separation, 
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Chh(,f)j is a test of clustering theory. Together these indicators provide valuable information as to 
whether galaxy groups have approximately the properties expected on the basis that the galaxies 
are orbiting in a dark-matter dominated potential that is at least partially virialized. 

Our statistical goals for group selection are to maximize the number of real groups and minimize 
the number of groups which contain rcdshift space interlopers. The derived group attributes come in 
three categories, with increasing sensitivity to the search parameters. First is the group centers, x, y 
and z; second are the group extensions in Vp, and velocity, and third there are the specific galaxies 
in the groups. For all purposes the locations are key, whereas the extensions are secondary. The 
extensions do enter scaling relations, and the precise group membership is not relevant for averaged, 
background subtracted measurements. 



4.1. Virial Mass-to-Light Ratio 

The ratio of the virial mass to the total luminosity, L^^, is a valuable indicator of the CDM 
mass density of the universe. It does not measure any component of the mass that clusters weakly. 
The virial mass-to-light ratio of groups has quite considerable scatter simply as a result of both small 
number statistics and orbital projection. On the basis of dynamical simulations Heisler, Tremaine 
&; Bahcall (1985) found the dispersion in virial mass estimates as a result of these fluctuations 
is nearly a factor of two above and below the true mass, but this was comparable to other mass 
estimators. We calculate the virial mass, Mvt, using the galaxies in the groups with < 1.5r2oo 
and Av < 3ai, following the CNOC methods (Carlberg et al. 1996). That is, 

Mvt = ^<^lRh, (2) 

where the virial radius is evaluated for the galaxies identified as being group members. The circu- 
larly averaged harmonic radius is 

\ i / i<j \ 3) 

where k^^ = 4Tirj/[{ri + Vj)'^ + s^] and K{k) is the complete elliptic integral of the first kind in 
Legendre's notation (Press et al. 1992). The softening, s = 2 arcsecond, eliminates the divergence 
for galaxies at the same radii from the group center (Carlberg et al. 1996). The luminosity, i^^, 
is k-correctcd, evolution compensated, and includes an extrapolation of the luminosity function to 
allow for galaxies below the redshift dependent absolute magnitude cutoff. The evolution is taken 
to be at a mean rate of one magnitude per unit redshift. 

If group galaxies are drawn from a universal luminosity function and the ratio of dark mass 
to luminous mass is a constant then the median Myr/L should be constant. The spread of the 
distribution of values can be used as an indicator of the statistical reliability of the group selection 
procedure. In Figure 3 we plot the median Myr/L against the fractional difference between the 



-8- 



first and third quartile Myr / L values. In the same system clusters have Mvt /L = 380± 70 M0 / L0 
(Carlberg et al. 1996), where we have removed the CNOCl correction for the mean flattening of 
clusters. The median Myt/L increases with both r™"^ and r"'"^'. Groups selected with r™"^ = 
0.5/i~^ Mpc and large r™"^ have huge median M/L values and a large spread between first and 
third quartile values. Smaller r'^^^ lead to an increase in the spread and decrease in the mean 
Mvt/L. The origin of this decrease in Myr/L with size is at least partially a result of the internal 
M/L gradient within groups that we discuss below. Overall, the r^"^ = 0.25/i~^ Mpc groups 
with r™"^ < 7/i^^ Mpc have the desirable property that both the spread and the median of the 
distribution do not change too much with r^°-^. 

The derived properties of the groups have very substantial uncertainties because of the small 
number of galaxies with velocities. The errors in the velocity dispersion and the resulting correlation 
with the derived virial mass-to-light ratios are illustrated in Figure 4 for the r™"^ = 0.25/i^^ Mpc 
and r^"^ = 5h~^ Mpc groups. The large errors axe the dominant source of the very strong correla- 
tion between the velocity dispersion and the derived M/L, as a consequence of Myr as 3G~^(Tfr„, 
which accurately predicts the slope of the correlation visible in the figure. The inset diagram re- 
stricts the sample to groups with at least six members. In this case there is some support for the 
indication that the the mass-to-light ratio rises with velocity dispersion. 



4.2. The Number Density- Velocity Dispersion Relation 

The number of groups as a function of their line-of-sight RMS velocity dispersion is given 
in Figure 5 for a range of group search parameters. The median velocity dispersions for r^"^ = 
Mpc are 192, 229, 256 and 266 for of 3, 5, 7, and lO^"^ Mpc, respectively. Below 
100 km s~^ is the regime of individual galaxies, which reduces the number of groups. A small effect is 
that the velocity precision of the survey is about lOOkms"^, which artificially reduces the numbers 
of low velocity dispersion groups. At this stage we recall that large values of the search length in 
the redshift direction tends to include enough outlier galaxies that a few groups arc promoted into 
the high ai tail of the distribution. Given that high velocity dispersion groups also tend to contain 
the most galaxies (for a constant M/L), large groups are the easiest to find. Thus low membership 
groups with high velocity dispersions are likely to have erroneously large ai values. 

The Press-Schechter (1974) theory works well to describe the abundance of halos in n- 
body experiments and a range of observational data, including clusters over a range of redshifts 
(Carlberg Morris Ycc & Ellingson 1997) and low redshift groups (Girardi &: Giuricin 2000). We 
can compare the Press-Schechter prediction for the number of groups to our observations. Because 
this requires absolute numbers we will use the redshift range of greatest completeness, roughly 0.2 
to 0.45. 

To calculate the expected density we follow the procedures outlined in Carlberg et al. (1997) 
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using Um = 0.2, J7a = and as = 1.0. We calculate the mass- velocity dispersion relation as 

Mi.5 = 8.6 X lOV? s^ km-^ M©, (4) 

where o"i is given in units of kms""*^ and the mass is the nominal value inside a virialized l.5h^^ Mpc 
sphere, as is appropriate for rich clusters, to which we want to normalize these predictions. The 
equivalent top-hat radius that contains this mass at the mean density is 



Rl 8.430°-2f/(3-f) 



M- 



1.5 



6.97 X 10i4j)^/i-i Mq 



l/(3-p) 



(1 + /j-i Mpc, (5) 



where p ~ 0.64 is the rate of increase of mass with radius (White Efstathiou & Frenk 1993; 
Carlberg Morris Yee & Ellingson 1997) and fi^ is the value of J^m at redshift z. We evaluate 
the Press-Schechter relation as 

n(M(ai))dM = ^'J'Jjyl^ exp [-S^A^)/2A^] dM, (6) 

where 5c{z) = 0.15{127rf/^n'^-°'^^^/D{z,n), (Navarro, Frenk k White 1996) and A{rL) is the top- 
hat fractional linear mass variance in spheres of radius calculated using a fitted CDM spectrum 
(Efstathiou, Bond &: White 1992). To determine the number density in bins of velocity disper- 
sion we simply integrate over the relevant range of masses. Note that the normalization we have 
used automatically means that our group number densities will match on to the CNOCl clusters 
(Carlberg Morris Yee & Ellingson 1997; Borgani, et al. 1999; Girardi & Giuricin 2000). 

The Press-Schechter predictions of number density for a median redshift of 0.36 are displayed in 
Figure 5. The subsample is contained in a volume of 1.8 x lO^hr^Mpc^ (or about 50% more for a flat 
cosmology). Below 100 kms"^ the sample is missing many halos for two reasons. First, individual 
galaxy halos make up the majority of the halos in this regime. The velocity dispersion of an 

elliptical is about equal to that of our median group (which in itself suggests an evolutionary 
connection). Second, because our velocity accuracy is about 100 kms~^ low velocity halos are 
scattered into the next higher bin. At this stage it we recall that our expected completeness rate 
for higher velocity dispersion groups is about 50% and even lower for those with velocity dispersions 
comparable to individual galaxies. 

The Press-Schcchtcr predictions are in reasonable agreement with the groups for r!^"^ = 
3h-^ Mpc and 5/i"^ Mpc, bearing in mind that the random errors are at least \/^. Smaller values of 
^max jj^jgg ]^[gii groups, while larger values, r™"^ > 7h~^ Mpc, produce a few highly improbable 
groups with the velocity dispersions of rich clusters. 

There are three reasons to select the r™"^ = Q.2hh^'^ Mpc and r™"^ = 5/i~^ Mpc groups as the 
best suited to our analysis of virialized halos with velocity dispersions of approximately 100-300 
kms~^, although other group selections parameters give rise to samples that show a very similar 
set of x,y,z locations. The selected groups have a relatively low dispersion in their M/L values, 
their n{ai) distribution is close to the Press-Schechter prediction with few high velocity outliers. 
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and the dynamical analysis below finds that our chosen redshift distance inclusion length pulls in 
most of the group members so that the derived velocity dispersions are fairly stable against the 
addition of more outlying members with increasing cutoff velocity. 

The locations on the sky of the "standard" r^"^ = 5h~^ Mpc and r^"^ = 0.25h~^ Mpc groups 
for one of the patches are shown as the points in Figure 6. The circles indicate the r2oo radii. 
Note that some groups are quite compact with respect to this radius. The parameters of all of the 
standard groups are given in Table 1. The columns give the location with co-ordinates measured 
relative to the designated group centers (Yee et al. 2000), ai (from which r2oo is calculated), the 
virial mass-to-light ratio, the number of group members with redshifts, and the mass. The random 
errors of the derived quantities, as estimated using the Jackknife technique, are very large for most 
of the groups, a straightforward consequence of the small numbers of members. 



4.3. The Two-point Group-Group Correlation Function 

A fundamental prediction of hierarchical dark matter clustering is that clustering, as measured 
by the two-point group correlation function ^GGif), should increase with the mass or velocity 
dispersion of the halo (Kaiser 1984; White et al. 1987). Here we measure both the redshift space 
correlation, ^(s),and the projected correlation function, Wp{rp), both of which provide an indication 
of the correlation length, tq. The co-moving redshift space separation is, 

with r(z) being the co-moving distance at redshift z. At separations small compared to the pairwise- 
velocity converted to a distance, ai2/H{z), ^(s) ceases to increase with decreasing separation as 
the random velocities begin to dominate the redshift space separation. On larger scales ^(s) is 
expected to be enhanced relative to ^ (r) as a result of the "compression effect" of systematic infall 
(Kaiser 1987). 

The projected correlation function, Wp{rp) = J ^[^Jvp + r^] dr^, has the advantage that 
the peculiar velocities have no effect on the result (Davis & Peebles 1983). We use the clas- 
sical DD/DR — 1 estimator (Peebles 1980) which is suitable in the strong correlation regime. 
The random distribution is derived from a fit to the observed n{z) distribution of the groups. 
Pairs are included for up to 30h~^ Mpc, co-moving. For a power law correlation function, 
^(r) = (ro/r)"', the reduced projected correlation function, Wp/rp is equal to A{'y)^{r), where 
A{j) = r(i)r(i(7 - l))/r(i7), a factor of 3.68 for 7 = 1.8 (Davis & Peebles 1983). 

We evaluate the correlations using the same procedures and programs used in Carlberg et 
al. (2000) to measure the correlation of galaxies. Using the standard group sample (r™"^ = 
0.25/1"-'^ Mpc and r^"^ = 5^~^ Mpc) we evaluate ^(s) and Wp{rp)/rp over the redshift range 0.15 
to 0.55. We use as much redshift range as possible to boost the sample size. The resulting 
correlation functions are displayed in Figure 7. In the form plotted both functions are dimensionless 
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functions. Note that the upward offset of Wp{rp)/rp relative to ^(s) is a natural result of the ^(7) 
factor. The error flags displayed in the figure are the square root of the number of galaxy pairs in 
each bin. Fitting the measured redshift space correlations to the function S,{r) = (so/r)^'*^, gives 
So = 6.8 zb 0.3/i~^ Mpc. The fit for the projected correlation function finds ro = 6.5 ± 0.3h~^ Mpc. 
These errors are formal fitting errors and do not include an allowance for the patch-to-patch variance 
which likely dominates the random error. 

The correlations measured for the CN0C2 galaxies over this redshift range in this cosmol- 
ogy have a mean of ro = 4.2 it 0.2^~^ Mpc where this error includes the patch-to-patch variance 
(Carlberg et al 2000b). The ratio of the correlation amplitude of our groups to that of the galaxies 
is 2.2 ± 0.4. At low redshift, Ramella, Geller & Huchra (1990) found sq ~ Sh'^ Mpc for the CfA 
groups, which are somewhat lower mean internal density than ours. A similar result emerged for 
the combination of the CfA and SSRS2 groups (Girardi Boschin &: da Costa 2000) which found 
that the group correlation amplitude was a factor of 1.64 ± 0.16 stronger than that of galaxies, 
although the CfA groups alone appeared to have a smaller offset. 

The biasing theory of Mo &; White (1996) gives predictions of the expected ratio of correlation 
amplitudes. The characteristic velocity dispersion of normal galaxies is about 100 kms~^, whereas 
our groups have a median velocity dispersion of about 200 kms~^. Prom the CDM power spectrum 
(Efstathiou, Bond & White 1992) with a shape parameter of F = 0.2, we calculate the linear mass 
variances to be 3.3 and 2.4, for 1 and 2h~^ Mpc top-hat perturbations, which are approximately 
the unperturbed mean radii of the perturbations associated with galaxies and groups, respectively. 
With Sc — 1.68 in the Mo &; White model, we find that groups should be more strongly correlated 
by a factor of about 1.75 or 1.47 with the Jing (1998) n-body calibrated modification. Provided 
that the group velocity dispersion is about twice the galaxy velocity dispersion the results are not 
very sensitive to the velocity dispersions chosen. We conclude that the correlation amplitudes are 
consistent with the expected relation at about the 1.5 standard deviation level. 

At this stage we have shown that our sample of groups is in good accord with three global 
statistics predicted on the basis of galaxies tracing cold dark matter halos. The next issue is to 
examine the relative internal distribution of the dark matter relative to the galaxies. 

5. Mean Internal Structure of the Groups 

N-body simulations now have sufficient length resolution that they can reliably predict 
the highly nonlinear realization of the collapse of halos and their resulting internal density 
profile. A fundamental prediction is that CDM halos have a central density cusp, roughly 
r-^ ° to (Dubinski & Carlberg 1991; Navarro, Prenk & White 1996; Moore et al. 1999a; 

Avila- Reese, Firmani, Klypin & Kravtsov 1999). At large radii the density begins to drop as ap- 
proximately r~^. The characteristic radius of the density profile can be calibrated in n-body 
experiments and derived from an approximate theory (Navarro, Frenk & White 1996). For the 
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massive dark matter halos of rich clusters the NFW density function, p{r) = Ar~^{r + a)~^, is 
entirely consistent with the derived mass distribution (Carlberg et al. 1997). The situation for 
galaxy mass halos is somewhat controversial, with constraining data coming from disk rotation 
curves in the presence of substantial amounts of baryonic mass and the modeling complication 
of partial pressure support relative to the circular velocity. However there is evidence that in 
small velocity dispersion halos CDM may allow central densities that exceed the observations 
(Moore et al. 1999a; Moore et al. 1999b). 

A fundamental difference between a galaxy and a rich cluster of galaxies is that galax- 
ies have a strongly rising mass-to-light ratios with increasing distance from the center whereas 
rich galaxy clusters have a nearly constant mass-to-light ratio over their virialized volume 
(Carlberg Yee &: Ellingson 1997). Therefore a basic question is whether groups exhibit a rising 
mass-to-light ratio. Individual groups have too few galaxies to make such a measurement. More- 
over groups come with a wide range of galaxy contents and have somewhat uncertain virialization 
because of their small numbers. Therefore we assemble a "stacked" mean group to boost the 
numbers to levels where we can make reliable measurements of the density and velocity dispersion 
profiles. 



Precisely how the mean group is built from the individual groups will have a significant impact 
on its properties. For instance, poorly determined centers will create a core in the projected profile, 
or, overlaying a small low velocity dispersion group with a large high velocity dispersion group will 
lead to a rising velocity dispersion with radius. The idea is simple: we stack galaxies in both Vp 
and Av on the group centers keeping track of the expected number of background galaxies. All 
galaxies are used, with no distinction given to galaxies that were used to define the group center. 
The group centers are determined in the virial analysis as the geometric weighted mean of the 
locations of the iteratively selected group members. We have adopted geometric weights, which 
help to compensate for the objects without redshifts. The smoothing radius used to calculate our 
geometric weights is 120" which is so large that there is little weight variation within these small 
groups, but weights do vary from group to group. That is, at the x,y,z location of each of the no 
groups in turn, we count the number of galaxies in the neighborhood at separations rp/r2oo Siiid 
Av/ai, to measure the stacked density distribution in redshift space of a group, ncgivp, Av). The 
stacked density distribution is related to the two-dimensional group-galaxy correlation function, 
^Gg as. 



where ng{z) is a smoothed fit to the redshift distribution locally converted to velocities in precisely 
the same way as the real galaxy pairs. This is the DD/DR—1 method (Peebles 1980) of calculating 
a cross-correlation between groups and galaxies. The correlation function ^Og is the fundamental 



5.1. The Stacked Mean Group 
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quantity of interest, since it has the background density removed. Both the surface density profile 

and the velocity dispersion profile are derived from ^Qg. Note that the total luminosity of the group 
uses the magnitude weights, which allows for the declining sampling rate towards the sample limit. 
Figure 6 builds the confidence that the selected centers are entirely sensible, although at this stage 
we cannot claim that this is an optimal procedure. 

To assemble our "standard" mean group we will scale all the velocities to the ai measured 
for each group and the radii to the r2oo derived from the velocity dispersion. Because there is a 
substantial uncertainty in ai for these small groups there is a concern that these scalings contain 
a large random element. At some level this is unavoidable. As the galaxies orbit in the groups 
they will have different velocities and positions with time which because of the little averaging 
available in small groups leads to considerable variation of their redshift space properties. Scaling 
the velocities and the radii avoids the pitfall of overlaying small, low velocity dispersion groups 
on the inside and large, high velocity dispersion groups on the outside, which would immediately 
lead to a rise of velocity dispersion with radius. The plot of the radial location of each galaxy that 
contributes to the mean group as a function of the source group's velocity dispersion in Figure 8 
does show one potential systematic problem: the higher velocity dispersion groups do not extend 
out to 1.5r2oo- This suggests that redshift space interlopers are an important source of noise leading 
to velocity dispersions that are too large for the size of the group. To address this issue we will be 
considering a number of alternate groups samples to check our results. 

Although these groups have small numbers of members interlopers can be accurately removed 
statistically. We use the smoothed ng{z) that we derive from the galaxy sample as a whole using 
the procedures of our clustering analysis (Carlberg et al 2000b). The random data are scaled in 
precisely the same way at the true data to determine whether or not they are in the cluster. To do 
the subtraction we need to bin the mean group in and Av bins. We select bins that arc 0.2 dex 
in logio(rp/r2oo) and 0.1 in the velocities scaled to ai. The binning is needed in our subsequent 
analysis anyway. We make velocity distance cuts at 3, 4 and 5 units as a check on the effect of 
residual interlopers and normally use the minimal noise 3 unit cut as our standard dataset. 

Although we have statistically subtracted the background, the velocity dispersion estimate is 
very sensitive to the presence of field interlopers. In physical co-ordinates the densities of the groups 
are expected to be about lOOpo near r2oo, but this is reduced by a factor of roughly 3ai/[H{z)r2oo] — 
20 at our At; = 3ai cutoff. Consequently interlopers, which axe distributed nearly uniformly 

in distance, generally cause erroneously large velocity dispersions. If there is a constant density 
normalization underestimate then the fraction of interlopers will increase with rp, since the contrast 
between the group and the surrounding field declines with increasing r^. 

The probability of group interlopers can be estimated from the two-point correlation function. 
The probability that a galaxy with Vp < r™'*^ and At; < Af™"^ is physically within r < rp is 
discussed in Carlberg et al. (2000a) . The conclusion is that slightly more than 50% of the galaxies 
within the required velocity separation are within the required radial separation, r. The other 50% 
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are clustered towards the group but most probably at roughly a correlation length. Other than 
small factors of order unity this redshift to real space ambiguity always exists. This means groups 
of size n may well be of size n — 1. The chances that they are of size n — 2 drops precipitously. 
Of course as n rises the fraction of interlopers is falling. These correlated but unvirialized group 
galaxies cannot be avoided. Although they present a challenge for a particular group they are not 
a significant complication for most statistical analyses of the mean group, given the mean group 
density profile we derive below. 



5.2. The Mean Density Profile 

Summing the mean group over velocities gives the mean projected density profile displayed 
in Figure 9. The points are for cutoffs at 3, 4 and 5 velocity units, with the points declining 
slightly at small radii for increasing velocity cutoffs as a result of the normalization of the total 
integrated density to unity. The errors are evaluated from the patch-to-patch variance. There is 
little systematic change with increasing velocity cutoff. It is immediately clear that the projected 
density profile is very nearly a power law with only a weak break to a slightly shallower central 
profile. The mean projected slope of our virialized groups is approximately Sjv oc R~^'^. The effect 
of restricting the sample to the ai < 200 km groups is shown in the inset to Figure 9. The 
density profile becomes somewhat steeper at large radii and there is a suggestion there is a break 
in the power law at small radii. The best fit parameters are a ~ 0.7 and b ~ 2.5, although the core 
slope is not well defined. 

We model the projected galaxy density distribution as the projection of the galaxy number 
density distribution, 

^(^) = 4nr-(r + cf' 
which is projected to a surface density using 

En{R} = 2 r u{r) dr. (10) 

so that we minimize the variance in the Y1n{R) plane. The resulting best fit is degenerate because 

the scale radius is always found to be very small. For the data displayed we find c = 0.061,0.074 
and 0.062, for cutoffs of 3, 4 and 5, respectively. The [a,b] pairs are [0.50,2.05], [0.81,1.71], and 
[0.78, 1.71], respectively. However, other fits with a + b = 2.55 ± 0.05 are equally acceptable given 
the small implied scale radius, c. 

It is somewhat remarkable that using our mean group center we identify such a cusped density 
profile. The procedure is not guaranteed to have any galaxy at the center of the group, whereas 
we find that on the average the galaxy density declines quite steeply away from the center. This 
is quite different from the cluster situation where a similar average center gives a reduced central 
density (Carlberg Yee &; Ellingson 1997). The NFW profile does not provide a good fit to this 
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distribution of galaxies, although it must be borne in mind that the NFW model describes the 
dark matter profile, not the galaxy distribution which can in principle be quite different. Our 
observed group profile is near a power law which is steeper in the center and shallower at large 
radii, than the NFW function or the steeper central cusps found by others (Moore et al. 1999a; 
Avila- Reese, Firmani, Klypin & Kravtsov 1999). 

5.3. The Mean Velocity Dispersion Profile 

The velocity dispersion profile allows the mass profile to be constrained. The projected velocity 
dispersion profile, ap{R), for the r^°^ = 5h~^ Mpc groups is shown in Figure 10. It is crucial that 
the weak rise, or at least absence of of a decline, in the projected velocity dispersion with radius 
be real, not an artifact of inadequate background subtraction. We display the results for several 
alternate values of r™"^ in Figure 11, and find that all give similar velocity dispersion profiles. 
Figure 11 also shows that setting r™**^ larger than the standard tends to increase the noise, while 
smaller values tend to leave out genuine group galaxies in the ai estimate so that larger velocity 
cutoffs quickly lead to higher ap{R). 

To help understand the distribution of possible outliers we present Figure 12, a gray scale of 
the Vp — Av plane of the mean group profile, normalized to ^Gg(^p) ^^)/^(''p)) fo remove the surface 
density variation. The plot shows that inside of about 0.05r2oo the statistics are too poor to give 
useful results. Beyond about 2r2oo there is no expectation that the galaxies are in virial equilibrium. 
The figure provides evidence that the slow rise in velocity dispersion across the virialized region 
is real and not a consequence of increasing interlopers with radius. Beyond about 2r2oo virialized 
motions are small and the width of the distribution increases rapidly. The slow rise of ap{R) is 
clearly rooted in the data, although a constant ap is not strongly excluded. It is well known that the 
pairwise velocity dispersion of galaxies has a similar weak rise with increasing separation; however 
that is expected on the basis of the two-point correlation function falling as r~^'^, rather than the 
r~^'^ that we find here. A dynamical constraint from the Jeans equation is that in a scale free 
distribution, oc p{r)r^. Therefore, the surprising outcome is the combination of the "steep" 
S(i?) and "slowly-rising" ap{R). Interlopers are always a concern but will tend to increase the 
outer values of both these functions, to leave the discrepancy in place. 

To evaluate the significance of the radial trends we require reliable errors estimates. The 
errors in S(-R) and ap{R) are robustly and empirically evaluated from the patch-to-patch variances 
in the quantities of interest. We find that the sensitivity to the velocity cutoff is relatively small, 
except that the errors increase in size. As a consequence we will use the 3 velocity unit cutoff as 
our standard. The velocity dispersion is calculated as the second moment of the binned projected 
velocity distribution function (pvdf). Figure 13 displays the average pvdf, normalized to the velocity 
dispersion in each radial bin. We note that the observed pvdf is much closer to Gaussian than to 
an exponential, although relative to the Gaussian the group mean pvdf appears to have a small 
excess at small and large velocities. This pvdf indicates that only mild anisotropy is allowed 



-16- 



(van der Marel et al. 2000), at least for simple distribution functions. 

The projected velocity dispersion is easily modeled as the projection of a 3D velocity dispersion, 
where we adjust the radial velocity dispersion, (Tr(r), and the velocity anisotropy, /3(r), to fit the 
data. The projection integral is, 

ctI{R)T.n{R) = v{r)al (^1 -/3(r)-^J ^ dr. (11) 

N-body models suggest that the orbits of galaxies should be similar to that of the full dark matter 
distribTitioii, with some deficiency of radial orbits which are selectively removed due to merging 
and tidal destruction (Ghigna, et al. 1998). In principle, the velocity ellipsoid is constrained by 
the pvdf, shown in Figure 13. 

We model the radial velocity dispersion as 

which is a solution of Jeans equation which worked well for clusters (Carlberg Yee & Ellingson 1997). 
In the absence of other information, the velocity anisotropy, /3 = 1 — o"g/o"^, is best taken to be a 
constant independent of radius. However, that choice immediately leads to the somewhat surprising 
result that the mass-to-light ratio rises with radius. In order to investigate the sensitivity of the 
result to the velocity anisotropy we adopt the model, 

/?(r) = /?oMr-) + /?oo[l-«^(r)], (13) 

where w(r) = r-|/(r| -|-r^). This /3(r) function goes to /3o at the origin and /3oo at large radius. As 
reasonable alternatives we set = 0.1r2oo or 0.3r2oo- This parameter could in principle be part of 
the nonlinear fit, but because dynamical problems are under-constrained without some knowledge 
of the velocity anisotropy we choose to use values of this parameter where our /3(r) gives a useful 
variation over the range of the data. The fits are shown for the three different velocity cutoffs in 
Figure 10. Note that the results for different /3 are essentially indistinguishable, and that the model 
has little difficulty giving a statistically acceptable fit to the data. 



6. The Rising Mass-to-Light Ratio 

The dynamical mass profile can be derived from the tracer density profile, v{t)^ its velocity 
dispersion, crr(r), and the velocity anisotropy, /3(r), using Jeans' equation. 



M{r) = 



d\a.a^ dlnv 



+ 



d In r d In r 



+ 2(5 



(14) 



The validity of this equation does not rely on v{r) being distributed like the mass density, p{r). 
It does require that the system be in dynamical equilibrium, which our group selection procedures 
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are specifically designed to pick out. Throughout we will use L(r) = 47r / ur'^ dr as the luminosity 

profile, assuming galaxies have the same mean luminosities at all radii. Strictly speaking we are 
working out the mass-to-number ratio. The assumption that low and high luminosity galaxies are 
similarly distributed relative to the mass field is explicitly tested below. The advantage of this 
approach is that it weights the galaxies relatively equally, rather than concentrating most of the 
statistical weight on the high luminosity galaxies. 

In Figure 14 we display the M{r)/L{r) profile for the standard group centers with the velocity 
dispersion calculation cutoff set at three velocity units. Note that this shows the integrated interior 
M/L, not the local values. We first model these standard groups using a conventional approach 
with minimal anisotropy. That is, the central velocity ellipsoid is isotropic, i.e. f3o = 0, the 
anisotropy radius, r^, is at about the midpoint of the virialized region, rp = 0.3r2oo- The outer 
anisotropy, P^o: takes on the values —1,-^,0 and ^, which range from nearly tangential to slightly 
radial velocity anisotropy. The M/L curves are normalized to unity at the data point closest to 
r/r2oo = 1- The data indicate about a factor of three to ten rise in the mass-to-light per decade 
of radius, with little sensitivity to the outer velocity anisotropy. That is, the light is much more 
centrally concentrated than the mass. This effect is visible for many of the groups where their sky 
distribution is much more concentrated than their r2oo as shown in Figure 6. 

The insets in Figure 14 address the concern about erroneously large ai values biasing the 
properties of the mean groups. The lower inset shows the results of the analysis based on the 
stacked groups restricted to have ai < 200kms~^. The upper inset is the result of the analysis 
where we have overlaid all groups in physical co-ordinates. For convenience and comparability in 
this case we have taken cri = 200kms~^ for the scaling. Although there are differences in the 
details of the result M/L curves the rising trend is preserved. Overlaying groups in physical space 
serves to minimize the size of the rise. The differences between these results serves as an indication 
of the size of the systematic errors. 

The mass-density profile, p{r), derived from the mass profile as (47rr^)~^dM/dr, is plotted in 
Figure 15 for the same models as in Figure 14. The model with the most extreme velocity ellipsoid 
at large radii, the tangentially dominant /3oo = — 1 (solid line), implies a large core which becomes 
relatively smaller as Poo becomes more positive. This model dependent result is one of the most 
interesting aspects of these groups and its reality needs to be tested which we can at least partially 
do within our own data. The upper and lower insets are for the ai < 200kms~^ sample and the 
physical co-ordinate stacking of the entire sample, respectively. 

An alternate version of the groups, with r^"^ = 3h~^ Mpc, is analyzed in the same way as the 
^max _ 5/j-i y[pc groups with mass-to-light shown in the upper left panel of Figure 16 and density 
in Figure 17. In this case, an even steeper rise of mass-to-light emerges than for the standard 
groups. For velocity ellipsoids with P^o > — 1 the central mass density is negative. The effect of 
analyzing the standard groups with rp = 0.1r2oo is shown in the upper right panels of Figures 16 
and 17. This small rp boosts the sensitivity of the gradient to Poo- We conclude that the existence 
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of a core in the dark matter distribution depends sensitively on the assumption that the velocity 
anisotropy is close to isotropic. 

A strong tangential anisotropy in the center, /3o = —1, is able to almost eliminate the mass- 
to-light gradient as shown in the lower panels of Fig 17, with = 0.3r2oo (left) and r/3 = 0.1r2oo 
(right). However we recall that the pvdf of Figure 13 does not favor such strong anisotropy. The 
implied density, shown in the lower left panel of Fig 17, is nearly a power law, r~^, for = — 1, and 
none of these resemble the predicted dark matter density profiles (Navarro, Prenk &: White 1996; 
Moore et al. 1999a). 

The source of the rising M/L in the kinematic data is readily understood. In a power law 
density distribution with a radial density dependence of p oc r~^+^, the mass generating the po- 
tential must have a velocity dispersion profile of cr^ oc r^. Prom the velocity dispersion we measure 
that p ~ I lb |. This implies a mass density profile that is slightly shallower than r~^. The 
measured mean light profile, z/(r) oc r^^^i^ with q ~ — ^. Consequently the mass to light ratio, 

3 1 

f pr'^dr/ J ur'^dr, which varies as r^^'^, rises as r4^4. Nearly circular orbits can flatten the M/L 
profile for the observed velocities, but these (3 would be very different than the observational results 
for rich galaxy clusters (Carlberg Yee & Ellingson 1997; van der Marel et al. 2000) and n-body ex- 
periments (Ghigna, et al. 1998). One significant difference between groups and rich clusters is that 
the timescale for dynamical friction to act declines as af, so that dynamical friction is much more 
effective for galaxies orbiting in groups than clusters. 



6.1. Color and Luminosity Dependent M/L Profiles 

As a test of the galaxy population dependence of the result, we divide the group galaxies in 
to blue and a red subsamples, splitting at (B — R)q = 1.25 mag. In a second test we divide the 
group galaxy sample into high and low luminosity subsamples, splitting at M^'^ = —20 mag. The 
resulting M/L profiles, both masses and luminosity profiles calculated from these subsamples alone, 
are shown in Figure 18. To simplify the comparison we keep (3o = for all models, but continue to 
vary /3oo- 

The blue galaxies do a much better job of tracing the mass profile, that is, it is fairly flat, 
than the red galaxies do. The opposite is true in clusters (Carlberg Yee Sz Ellingson 1997). In 

both clusters and groups the red galaxies arc more centrally concentrated than the blue galaxies. 
The M/L gradient does not have a signiflcant dependence on the luminosities of the galaxies. The 
latter is an argTimcnt that dynamical friction on a galaxy's gas and stars might not be the dominant 
source of the segregation of mass and light in groups. However, the presence of a much stronger 
M/ L gradient for the red galaxies would be consistent with dynamical friction taking several orbits 
after entry into the cluster to produce a segregation. 
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7. Discussion and Conclusions 

In this paper we use the CN0C2 survey catalogues, containing approximately 6000 galaxies 
with redshifts, to identify more than 200 high probability virialized groups containing, on the 
average, 3.8 members with redshifts for a total of about 1000 group members. This represents 
about 25% of all eligible galaxies in the sample. These groups have selection dependent Myr/L 
values comparable to the values of large clusters with the "standard" groups being about a factor 
of two lower. The group population volume density matches on to the number-velocity dispersion 
relation of clusters. The clustering of groups is enhanced at about the level expected for the 
approximately 200 kms~^ velocity dispersion dark matter halos of groups as compared to the 
approximately 100 kms~^ halos of individual galaxies. Overall the global properties of groups are 
about as expected for dark matter halos on this mass scale. 

Quite unlike rich clusters which have a mass-to-light profile constant with radius 
(Carlberg Yee & Ellingson 1997), our analysis of the internal properties of galaxy groups finds 
considerable evidence for a gently rising mass-to- light profile with radius. A second difference is 
that in groups the blue galaxies follow the dark mass more closely than the red galaxies. This paper 
presents several lines of argument that the rising velocity dispersion profile with radius is not an 
artifact of the analysis. Moreover, the projected velocity distribution function indicates that the 
orbital distribution is unlikely to consist almost entirely of circular orbits, as would be required to 
minimize the rise in M/L with radius. 

The inference that M/L rises with radius depends sensitively on the theoretical and indirect 
observational indication that the mean velocity ellipsoid is approximately isotropic. Hence this 

result needs to be independently verified. Current indirect support for this result comes from X- 
ray observations of low redshift galaxy groups having velocity dispersion overlap with our sample 
near 300 kms~^. The X-ray emission is significantly more extended than most of the group 
light indicating a rising M/L (Kriss Cioffi & Canizares 1983; Mulchaey & Zabludoff 1998). Weak 
gravitational lensing is one of the best prospects to check this result since it has no dependence 
on assumptions of equilibrium. These intermediate redshift groups are ideally situated for weak 
lensing studies (Hoekstra et al. 2000). The somewhat puzzling weak lensing results for the low 
velocity dispersion cluster MS1224-I-20 in two independent studies indicates that the M/L is rising 
with radius (Fahlman Kaiser Squires & Woods 1994; Carlberg 1994; Fischer 1999). Given the large 
group to group fluctuations it will be important to examine a statistically meaningful sample. 

A highly model dependent result is that we find that the mean group most likely has a signif- 
icant core radius. To retain a power law form would require very tangential orbits near the center. 
However, a large core could simply be a trivial result of the group galaxies sloshing around in the 
group potential, such that the centers of light and mass are not coincident. To further investigate 
the location and radial dependence of the central density profile of groups will require studying 
individual groups in detail, possibly with X-ray and weak lensing measurements. 

What is the origin of the rising M/Ll The two physical possibilities are that galaxies sink 
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with respect to the dark matter, or, that group dark matter has a small effective pressure which 
does not allow the core to undergo the same gravitational collapse as cold collisionless matter can. 
It is important to note that near r2oo both the galaxies and the dark matter have undergone about 
the same amount of collapse relative to the field. It is only in the inner third or so where the 
major differences develop, and the size of those differences is strongly dependent on the kinematic 
model. The simplest explanation is that these observations confirm the standard picture of a 
collisionless dark matter with gaseous dissipation allowing the baryonic matter to concentrate to 
the center of individual dark matter halos. As galaxies cluster together to form a group their 
individual halos are tidally removed to join the common halo from the outside in. Dynamical 
friction from the collisionless dark matter causes the galaxies to sink in the common group halo, 
over about a Hubble time (Barnes 1985; Mamon 1987; Evrard 1987; Bode Cohn & Lugger 1993; 
Pildis Evrard &; Bregman 1996). One possible concern with this interpretation is that when groups 
join together to make rich clusters this history needs to be erased to leave no measurable M/L 
gradient. Furthermore, whether this mechanism is consistent with no luminosity dependence of the 
M/L gradient but a large galaxy color dependence of the M/L gradient puts some fairly strong, 
but not necessarily unreasonable constraints, on the formation history of the galaxies. 

A more extreme possibility for the origin of the M/L gradient is that the dark matter 
is subject to some effective pressure that does not allow it to undergo full gravitational col- 
lapse to form a core, either through a phase density limit, or through collisional interactions 
(Spergel &; Steinhardt 1999). More conventional explanations should be carefully examined and 
our results independently verified with other methods of observation before this is accepted. 
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Fig. 1. — The line of sight velocity dispersion as estimated on the first pass (abscissa) of the group 
refinement and on the fourth pass (ordinate) for sets of group finding parameters. The r^°^ is 
fixed at 0.25/i~^ Mpc. The asterisks, circles, x's and squares are for r™"^ of 3, 5, 7 and 10h~^ Mpc, 
respectively. 

Fig. 2. — The rcdshift distribution of groups found with r™"^ fixed at 0.25h~^ Mpc and r^"^ 
(indicated on the figure in the upper right) varying over the same range as Figure 1 . 

Fig. 3. — The median M/L value of the groups versus the factional spread between the first 

and third quartilc M/L values. Squares indicate r™"^ = 0.5h~^ Mpc groups and circles indicate 
^max _ Q 25/i-i Mpc groups. The numbers in the symbols give the values of r^"^. 

Fig. 4. — The mass-to-light ratios of the groups as a function of their velocity dispersions for 
^max _ o.25/i~^ Mpc and r^"^ = 5h~^ Mpc. The la error flags plotted are Jackknife technique 
estimates. The correlation is the result of the af dependence of Myx- The inset restricts the 
sample to groups with six or more members. 

Fig. 5. — The distribution of measured line-of-sight velocity dispersions in lOOkms"^ bins, plotted 
at the bin centers. The parameters for the group selection are indicated by the symbols in the 
diagram. The dotted line shows the Press-Schechter prediction for the distribution using the cluster 
normalization. 

Fig. 6. — The xy locations of field galaxies (open circles) and group galaxies (filled), along with 
r2oo radii of the groups found between redshifts 0.1 and 0.45 in the 0223+00 field. The starting 
parameters are r™"^ = 0.25h~^ Mpc and r^"^ = 5h~^ Mpc (co-moving). The group members are 
marked with filled symbols. Note that the group members fill their r2oo to varying degrees. None of 
the members are outside the estimated region of virialization estimated for a low density universe, 
l-5r2oo- The major ticks are at intervals of 1000 arcsec, with North up and East to the left. 

Fig. 7. — The redshift space (open circles) and projected real space (filled diamonds) correlations for 
the r^"^ = 0.25h~^ Mpc and r^"^ = 5h~^ Mpc groups. The errors here are simply the square root 
of the pair count in the bin. The fitted correlation lengths are 6.8±0.3^~^ Mpc and 6.5ib0.3^~^ Mpc, 
respectively. 

Fig. 8. — The scaled radius versus the source group's velocity dispersion for all galaxies in the mean 
group for Au™"^ = 3ai prior to statistical background subtraction. The small plus signs are for all 
galaxies contributing to the mean group whereas the filled octagons mark the galaxies identified as 
group members using our algorithm. 

Fig. 9. — The number weighted surface density of galaxies in the mean group. The points and 
three fitted lines show the slow decrease of the normalized surface density at intermediate radii 
with clipping at 3, 4 and 5 velocity units. The inset shows 'En{R) using only the ai < 200 km s^^ 
groups with (solid line) and, for reference, all velocity dispersion groups cutoff again (dotted line). 



-30- 



both with 3 velocity unit clipping. 

Fig. 10. — The number-weighted projected velocity dispersion profile. The circle, cross and square 
are for velocity cuts at 3, 4 and 5 velocity unit cutoff, which are fit by the solid line and the dashed 
and dotted lines (which are hard to distinguish because they are overlaid). 

Fig. 11. — The number-weighted projected velocity dispersion profile for alternate values of r^"^ 
of 2, 3, 7 and 10 h"^ Mpc. The circle, cross and square are for velocity cuts at 3, 4 and 5 velocity 
units, and fitted by solid, dashed and dotted lines, respectively. 

Fig. 12. — The mean radially normalized group in position and velocity, ^(i?, At>)/S(rp), selected 
using r^"^ = 0.25h~^ Mpc and r^"^ = 5h~^ Mpc. Each column of fixed R is normalized to a total 
of unity. The gray scale is proportional to the logarithm of the local over-density in the normalized 
velocity and radius. The second moment of this distribution in the vertical direction is (Jp{R) of 
Figure 10. The lowest over-density value (white) is 0.2% of the total, and the highest (black) is 
10% of the total in the entire R bin. 

Fig. 13. — The projected velocity distribution function (pvdf), normalized to have unit velocity 
dispersion and unit area under the curve. The pvdf for cutoffs of 3 and 4 velocity units are shown 
as solid and short dashed stepped lines, respectively. The reference lines are a Gaussian and an 
exponential distribution, solid and short-dashed lines, respectively. 

Fig. 14. — The rising mass-to-light ratio from the 3 velocity unit cutoff velocity dispersion profile for 
group centers identified with r^"^ = 0.25/i~^ Mpc and r"''"^ = 5/i~^ Mpc. The galaxy co-ordinates 
are scaled to their groups ai and r2oo- The errors are computed from the velocity dispersion errors. 
The upper inset shows the results obtained by overlaying the galaxies in physical co-ordinates 
assuming a\ = 200kms~^ and the lower shows the results after restricting the analysis to groups 
with velocity dispersion less than 200 kms~^. The solid, dashed, dot-dashed and dotted lines are 
for /3oo = —1,-5,0, and |, respectively 

Fig. 15. — The inferred dark matter density profile for the variable velocity anisotropy model, with 

(5q = 0, rg = 0.3r2oo and /?oo = —1,-^,0 and \ for solid, dashed, dot-dashed and dotted lines, 
respectively, as in Figure 14. The lower inset shows the results obtained by overlaying the galaxies 
in physical co-ordinates and the upper shows the results after restricting the analysis to groups 
with velocity dispersion less than 200 kms~^. 

Fig. 16. — Profiles of M/L for an alternate data set, r™"^ = 3/i~^ Mpc and r™"^ = 0.25/i^^ Mpc 
(upper left) and the standard groups sample analyzed with a small rp = 0.1r2oo (upper right). In 
both lower panels we use an extremely tangential velocity ellipsoid, (3q = —1 with the standard 
rp = 0.3r2oo (lower left), and a small rp = 0.1r2oo, (lower right). The lines have the same meaning 
as in Figure 15. When the M/L has a zero-crossing dip the the implied values of the mass are 
negative at small radii. These show that a nearly tangential velocity ellipsoid can eliminate the 
inference of a rising mass-to-light ratio. 
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Fig. 17. — The density profiles for the data and models of Figure 16. Again, masses are negative 
at small radii when a zero-crossing dip is present. 

Fig. 18. — The inferred mass-to-light ratios (normalized to unity at r2oo) for galaxy subsamples 
within the groups. The top row has blue (left) and red (right) selected galaxies. The bottom row 
has high (left) and low (right) luminosity galaxies. The smooth lines show the modeled results. 
The main result is that blue galaxies trace the mass distribution much more closely than the red 
galaxies. 
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